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16. Abstract 

This report presents the results of a study in which the applicability at subsonic speeds 
of several aerodynamic methods for predicting dynamic gust loads on aircraft, including 
active control systems, was examined and compared. These aerodynamic methods varied 
from steady state to an advanced unsteady aerodynamic formulation. 

Brief descriptions of the structural and aerodynamic representations and of the motion 
and load equations ai*e presented. Comparisons of numerical results achieved using the 
various aerodynamic methods are shown in detail. From these results, aerodynamic 
representations suitable for dynamic gust analyses are identified. 

It was concluded that several aerodynamic methods are satisfactory for dynamic gust 
analyses of coafiguration.s having either controls fixed or active control systems that 
primarily affect the low frequency rigid body aircraft response. If an active control 
system is included that primarily affects the elastic modes response, the more exact 
unsteady aerodynamic methods are required to provide the necessary accuracy for 
predicting both the magnitudes and phase angles of the modal responses. 
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COMPARISONS OF SEVERAL AERODYNAMIC METHODS 
FOR APPLICATION TO DYNAMIC LOADS ANALYSES 


by Richard L. Kroll and Ronald D, Miller 
Boeing Commercial Airplane Company 


1.0 SUMMARY 


This report presents the results of a study in which the applicability at subsonic speeds 
of several aerodynamic methods for predicting d3mamic gust loads on aircraft, which 
include active control systems, was examined and compared. These aerodynamic 
methods varied from steady state to an advanced unsteady aerodynamic formulation. 
The study was performed under NASA contract NAS2-7729, ‘‘Development of a 
FLEXSTAB Computer Program,” (task IV, item 3). 

Brief descriptions of the structural and aerodynamic representations and of the motion 
and load equations are presented. Comparisons of numerical results achieved using the 
various aerodynamic methods are shown in detail. From these results, aerodynamic 
representations suitable for dynamic gust analyses are identified. 

It was concluded that several aerodynamic methods are satisfactory for dynamic gust 
analyses of configurations having either controls fixed or active control systems that 
primarily affect the low frequency rigid body aircraft response. These include the 
ILEXSTAB steady state aerodynamics'’’ modified with Kiissner and Wagner indicial 
functions, strip theory aerodynamics’'''’' modified with Kiissner and Wagner indicial 
functions, the Doublet Lattice unsteady method,*'" the vortex spline-kernel function 
unste'iii; method developed under this contract,"'* and the Rho unsteady kernel function 
aeroaynamic method.** The approximations to the more exact unsteady representation, 
that are pr . ided by the FLEXSTAB steady state and strip theory aerodynamic methods 
modified with Kiissner and '' lagner indicial functions, are suitable only because of the 
gust power attenuaLion at th^ hi;rh >r frequencies. If an active control system is included 
that primarily affects tbi , istie modes rc.sponse, only the Doublet Lattice, vortex 
spline-kernel funetinn, anu : ' itho unsteady kernel function aerodynamic methods are 
suitable, because ndy the . t t .’;act unsteady aerodynamic methods can provide the 
required accuracy in predi. LJug both the magnitudes and phase angles of the modal 
responses. 

The only methods currently capable of providing a complete aerodynamic representation 
for three-dimensional aircraft configurations are the FLEXSTAB solution modified with 
Kiissner and Wagner functions and the Doublet Lattice solution. The other aerodynamic 
methods are restricted to either planar surfaces and/or two-dimensional configurations, 
without including fuselage aerodynamics or interference effects. 


"Applicable to both subsonic and supersonic flow. 
’""Subsonic flow. 


2.0 INTRODUCTION 


In aircraft design^ differences in the unsteady aerodynamic methods used for predicting 
the interactive coupling between the aerodynamics and structural motions of lifting 
surfaces can dramatically affect the calculated stability of the system and, therefore, the 
calculated elastic response and structural loads. These methods must accurately predict 
both phase angle and magnitude of the pressures for each of the elastic modes. This 
requirement becomes less important at higher frequencies for dynamic gust analyses if 
controls are fixed or if the control motions are used only to influence the rigid body 
motions, because of the large attenuation of the gust pov^er at the higher frequencies. 
However, this requirement becomes critical with the inclusion of active control systems 
that are designed to alter the elastic structural motions in a presc ibed manner as in a 
modal suppression system. 

Every unsteady aerodynamic method is based on a numerical scheme containing 
approximations. The influence of these approximations on the predictions of dynamic 
loads and response of an aircraft is not adequately ascertained by simply examining the 
real and imaginary parts of predicted pressure distributions. Small differences in 
predicted pressure magnitudes and phase angles may dramatically affect the dynamic 
response and loads due to coupling effects between modes. To reliably assess the 
relative accuracy of particular methods in obtaining the interactions and resulting 
response, it is necessary to perform a full dynamic analysis of an aircraft including 
aerodynamic, structural, and control system representations. Comparisons can then be 
made with results predicted by previously validated aerodynamic methods. 

The purpose of this study was to compare the accuracy of several unsteady aerodynamic 
methods by using each method in performing dynamic loads analyses. The aerodynamic 
methods compared were: (1) an "s” plane'^' (Kussner-Wagner) modification of the steady 
flow method employed in the present PLEXSTAB program (ref. l); (2) the steady flow 
method plus the low frequency, constant pressure panel method of the existing 
PLEXSTAB program; (3) vortex spline-kernel function method developed under 
contract NAS2-7729, task IV, item 2 (ref. 2); and (4) the Rho unsteady kernel function 
method (ref. 3). In addition, results from a separate study (conducted under the Boeing 
IR&D program) involving two other aerodynamic methods are included for 
comparison.'*'*'^' They are the lifting line (strip theory) aerodynamics'*'*'*' modified with 
Kiissner and Wagner functions and the Doublet Lattice method of reference 4. The 
validity of all of these methods was measured using the Rho unsteady kernel function 
method (ref. 3), a current advanced state-of-the-art unsteady aerodynamic 
representation, as a basis for comparison. 


"^LaPlace transform parameter. 

‘^“^■Richard, M. and Kroll, R. 1., "Limitations of the Doublet Lattice Method-Application to a 
High Aspect Ratio Wing,” Boeing document D6-42395, 1975. 

'•*'^"=*Richmond, L. D., "A Rational Method of Obtaining Three-Dimensional Unsteady Aerodynamic 
Derivatives ofintersecting Airfoils in Subsonic Flow.” Boeing document D6-7401, 1962. 


3.0 SYMBOLS AND ABBREVIATIONS 


A 

AIC 

Ain> Aflut 

AR 

Afooti Atip 


ail ^>ii Cii di 
a 2 , bo, C2, d2 
br 


c.g. 

{Cl}, {C2} 

{C1UC2} 


P 

Fc 

fn 

g 

I 

k 

L 

M 

MAC 

ni 

Min, ^out> ^root 


gust response factor (o'l./o'g), response/unit velocity 
aerodynamic influence coefficient 

acceleration (g's), inboard and outboard nacelles, respectively 
aspect ratio 

acceleration (g’s), wing root and wingtip, respectively 

coefficients of Wagner function 

coefficients of Kiissner function 

wing root semichord = 326 in, 

center of gravity 

real and complex generalized forcing function coefficients, 
respectively 

real and complex coefficients of loads due to the forcing 
function, respectively 

forcing function 

compressibility factor 

natural frequency, cps 

acceleration due to gravity 

mass moment of inertia 

reduced frequency, , rad/s 

scale of turbulence, ft 

Mach number 

mean aerodynamic chord 

mass 

bending moment (in-lb), inboard nacelle, outboard nacelle, 
and wing root, respectively 


[M,] 

[Ma] 

[M4] 

[M5] 

[Ma] 

[M4] 

[Mg] 


No 

{P} 

[Prc]* [Plm] 
PSD 

{Q} 

q 

[Qrc]> [Qlm] 

rms 

T(i£i) 

V 

Vin> Vout, Vroot 


W 

w 

WBL 
X, y 


“B 


generalized structural stiffness matrix 

generalized structural inertia matrix 

generalized aerodynamic stiffness matrix 

generalized aerodynamic damping matrix 

inertia force coefficients 

aerodynamic stiffness force coefficients 

aerodynamic damping force coefficients 

number of zero crossings with positive slope/ft 

complex surface pressure matrix, force/unit area 

real and imaginary pressure matrix coefficients, respectively 

power spectral density 

complex generalized force matrix 

generalized coordinates 

real and imaginary coefficient matrices of the generalized 
forces, respectively 

root mean square 

complex load frequency response function 
airplane forward true speed, length/s 

shear (lb), inboard nacelle, outboard nacelle, and wing root, 
respectively 

downwash (w/v) 

velocity normal to lifting surface, length/s 

wing buttock line 

surface coordinates, length 

function of gust angle 



coefficients of Wagner function 

“2. ^2. 72 

coefficients of Kiissner function 

Ap 

surface pressure, force/unit area 

s 

damping ratio (c/cc) 

17 . ^ 

surface integration variables, length 

p 

density of air, mass/unit volume 

o-g 

rms of gust velocity, length/s 

o-r 

rms of response 

$(i£l) 

response quantity spectrum 

$(X1) 

gust spectrum 

*(t) 

Wagner function 


mode shape 

•Wt) 

Kiissner function 

ft 

spatial function (<u/v), rad/unit length 

0} 

frequency, rad/s 


iudicial convolution 

Subscripts 


Bdg 

bending 

S 

gust 

Im 

imaginary 

r 

response 

Re 

real 

Tors 

torsion 
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4.0 ANALYSIS 


In order to effectively compare each aerodynamic method employed in the dynamic 
analysis, the differences in the analysis results must be reflected only from the various 
aerodynamics. Thus, the structural model, structural modal interpolation routines, 
formulation of the equations of motion and load equations, and solution routines should 
remain fixed while only the aerodynamic methods are varied. 


Since the formulation and solution routines used in dynamic gust analyses involving 
active control system effects are based on linear analysis techniques, it is possible to 
separate the solution into two separate analyses: (1) dynamic gust loads analysis 
performed in the frequency domain with controls fixed and (2) dynamic loads analysis 
performed in the frequency domain with the excitation function consisting of an 
oscillatory control surface. The different aerod 3 mamic methods were employed in each 
analysis where applicable, thus, enabling comparisons of the ability of each 
aerodynamic method to model two basic types of dynamic analyses. The airplane 
configuration chosen for this study was a simplified model of a large subsonic jet 
aircraft. Other parameters used in relation to the study were: 

Speed 359 kn 

Mach number 0.8 

Altitude 20 000 ft 

Weight 753 000 1b 

C.G. location 16.2':?^^ MAC (MAC = 327.75 in.) 

Aspect ratio (AR) 6.96 

4.1 STRUCTURAL REPRESENTATION 

In order to simplify the structural model but maintain a meaningful dynamic analysis, 
the airplane structural model was represented by using only the wing elastic degrees of 
freedom. The body, empennage, and nacelles were considered to be rigid. The wing was 
modeled as a finite number of lumped masses, rotary inertias, and beam segments. The 
generalized mass and stiffness matrices were obtained for the first eight free-free modes 
of this model which consisted of two rigid body modes (vertical translation and pitch) 
and six wing elastic modes. This representation of the generalized mass and stiffress 
was used with each aerodynamic method in formulating the equations of motion. A 
more detailed description of the structural model is presented in appendix A. 

4.2 AERODYNAMIC REPRESENTATION 

The aerodynamic representation was dependent on the aerodynamic method employed. 
Following is a brief description of each method. 





4.2.1 FLEXSTAB AERODYNAMICS 


The panel-type aerodynamic method (ref. 1) in FLEXSTAB calculates either subsonic or 
supersonic steady state aerodynamic influence coefficient (AIC) matrices relating 
surface pressure to flow incidence, and a low frequency approximation to an unsteady 
AIC matrix relating surface pressure to flow incidence time rate of change. The 
aerodynamic theory considers linear potential flow in both subsonic and supersonic 
regimes. Solutions are written as integral equations involving the strengths of flow 
singularities distributed over aerodynamic mean surfaces and mean lines. Strengths of 
the flow singularity distributions are determined using approximations based on those 
used in the finite element method of Woodward. Thin body and interference body mean 
surfaces are represented by panels, each panel having a constant strength vortex flow 
singularity. The aerodynamic grid used for this method is shown in figure 1. 

The low frequency approximation to unsteady aerodynamics used in FLEXSTAB has 
severe frequency limitations and is generally suitable only for calculation of dynamic 
stability derivatives. The restrictions are based on reduced fre quency and an AR-Mach 
number relationship (k«l, k«(l-M^)/M^, k| Inkj <<2/AR yl-M^ for subsonic, and 
k«l, k«(M^-l)/M^ for supersonic flow). 

4.2.2 VORTEX SPLINE-KERNEL FUNCTION UNSTEADY AERODYNAMICS 

This panel-type method (ref. 2) is an unsteady adaptation of the steady vortex 
spline-kernel function method of reference 5. The integral equation relating pressure 
and normal wash on a lifting surface is solved using localized spline-type functions to 
represent the pressure distribution. These pressure functions are defined over a local 
portion of the aerodynamic surface, and a spline technique is employed to ensure 
continuity across adjacent panels. A modified version of the Rho unsteady aerodyuEimic 
method (ref, 3) for evaluating the kernel function is employed using these pressure 
functions with the spline technique. This finite element method is in its early 
development stages and is presently restricted to single planar surfaces in subsonic 
flow. The paneling scheme used for this method is shown in figure 2. 

4.2.3 RHO UNSTEADY KERNEL FUNCTION AERODYNAMIC METHOD 

The Rho unsteady kernel function aerodynamic method (ref. 3) solves the classic 
boundary value problem by equating the kinematical downwash to the integral 
downwash expression as shown in equation (1). 

ff^P(?.b) K(x, ?,y,i?,k.M)d|d77 (I) 

' 4jrpV^ 

First, the singularities in the integral expression are isolated and evaluated 
individually. Next, pressure loading functions that account for discontinuities in the 
downwash due to the control surface are assumed. These pressure functions are then 
substituted into the integral downwash expression with the proper discontinuities to 
calculate a mathematical downwash which is used to remove the discontinxiities in the 
kinematic downwash, resulting in a modified kinematical downwash that contains no 
discontinuities. A second set of pressure functions are then assumed which contain 
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unknown coefficients. These coefficients an* determined using the modified kinematical 
downwash in the solution. The total unsteady pressure loading becomes the sum of all 
the pressure functions. 

Although this method is presently the most exact theoretical formulation of subsonic 
unsteady aerodynamics available, it is restricted to single planar surfaces. This severely 
restricts its use in dynamic loads analysis where the total airplane must be represented. 

The Kho unsteady method is a collocation point method. The distribution of the 
collocation points used over the wing surface is shown in figure 3 , 

4.2.4 DOUBLET LATTICE 

The Doublet Lattice method tref, 4) evaluates the integral equations relating pressure 
and normal wash on lifting surfaces. Landahl's representation of a nonplanar kernel 
function is used. Streamlined closed-body aerodynamics are represented with Mile’s 
slender body theory. Woodward’s method of interference surfaces is applied to determine 
the lifting-surface slender-body interaction effects. The Doublet Lattice method is a 
subsonic panel-type method which can be used to represent the aerodynamics on the 
total airplane. The primary restiictions are that the aspect ratio of each panel should be 
approximately unity and that the number of panel divisions in the chordwise direction 
must increase with reduced frequency. Figure 4 shmvs the panel distribution used for 
the gust excitation problem; figure 5 illustrates the paneling scheme used in the 
oscillatory aileron case, 

4.2.5 STRIP THEORY (LIFTING LINE) 

This aerodynamic method^' is based on two-dimensional incompressible flow in 
formulating the aerodynamic lift curve slopes. These two-dimensional lift curve slopes 
can be modified with compressibility^ factors (F^J and wind tunnel test data where 
applicable. Steady state three-dimensional lift is obtained by modifying the 
two-dimensional lift coefficients with a static induction matrix based on the method 
developed in reference 6 to include three-dimensional effects. Strip theory aerodynamics 
are restricted to high AR lifting surfaces. The paneling used for this aerodynamic 
method is shown in figure 6. 

4.2.6 KUSSNER AND WAGNER LIFT GROWTH FUNCTIONS 
AS APPLIED TO THE STEADY STATE AERODYNAMIC 
FORMULATIONS OF FLEXSTAB AND STRIP THEORY 

One method of formulating the unsteady aerodynamic problem is to use the FLEXSTAB 
or strip theory steady state aerodynamic formulations with the Kussner and Wagner lift 
growth functions. Thi*- approach is discussed extensively in reference 7^, Briefly, the 
circulatory portion of the aerodynamic forces is obtained by multiplying the steady state 
lift and moment by the convolution of the normal wash (the velocity normal to the 


wf 

^Richmond. L. D., "A Rational Method of Obtaining Three-Dimensional Unsteady Aerodynamic 
Derivatives of Intersecting Airfoils in Subsonic Flow.” Boeing document D6-7401, 1962, 


lifting surface) with either the Kussner or Wagner lift growth function. Normal wash 
caused by structural motion is convoluted with the Wagner lift growth function which 
was derived assuming a step change in angle of attack. This convolution multiplied by 
the steady state aerodynamics gives the circulatory portion of the response 
aerodynamics. The cii’culatory part of the gust aerodynamic forces is calculated using 
the normal wash, resulting from a gust encounter, convoluted with the Kussner lift 
growth function. The Kussner function was derived assumming a sharp-euge gust. The 
exact form of the Kiissner-Wagner functions is not expressible in terms of well known 
functions; therefore, approximate representations are used; 


-Ofl 


0(t) = a j “ b| e 

-a-(t 




*’ -djc (Wagner) 

( 2 ) 

•yot 

-d-je ^ (Kussner) 

( 3 ) 


The values of the coefficients used in equations (2) and (.3) are dependent on aspect ratio 
and Mach number. References 8 and 9 present the development of the approximations 
and the tabular data for selecting the coefficient values. For this study, the coefficients 
were selected based on an aspect ratio of six and a Mach number equal to zero. The 
numerical values used for the Kussner and Wagner functions are shown in appendix B. 


4.3 EQUATIONS OP MOTION AND LOAD EQUATIONS FORMULATION 


The form of the equations of motion and load equations are presented in appendix B. In 
summary, two basic forms were used depending on the aerodynamic representation; 
(1) constant coefficient equations which are formed by utilizing steady state 
aerodynamic representations, and approximating unsteady aerodynamics with either 
Kussner and Wagner indicial lift growth functions or the low frequency representation 
of PLEXSTAB; and (2) nonconstant coefficient equations with respect to frequency, that 
are formed utilizing the more exact unsteady aerodynamic methods. The constant 
coefficient generalized structural inertia and stiffness forces were obtained from the 
structural representation and used in conjunction with both of the two basic forms of 
the aerodynamic force equation", The structural damping force was set equal to zero. 


Since it was the purpose of this study to obtain a comparison between several 
aerodynamic methods when used with a dynamical system, it was felt that the omission 
of structural damping would yield the clearest comparison of differences between the 
several methods in predicting both magnitudes and phase angle relationships. 
Admittedly, the addition of structural damping would have effect on the overall 
response of the system and, therefore, could change the magnitude of the relaiive 
differences between several aerodynamic methods calculated without structural 
damping. 

The load equation formulation used the force summation technique (ref. 10). Shears and 
bending moments were obtained at three wing locations; wing buttock line (WBL) 105, 
470, and 837, which correspond to wing span stations at approximately the wing root, 
inboard nacelle, and outboard nacelle, respectively. Accelerations were obtained at four 
wing stations; WBL 1132, 837, 470, and 64, which approximately correspond to wingtip, 
outboard nacelle, inboard nacelle, and wing root, respectively (fig. 7). 
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The acceleration equation coefficients are only a function of the mode shape. Shear and 
bending moments are a function both of the inertia load and aerodynamic forces. The 
inertia load coefficients were calculated using the mass and geometry data and then 
combined with the aerodynamic load coefficients to produce a total load representation. 

The methods of calculating the aerodynamic load coefficients for the different 
aerodynamic theories were similar with only slight differences caused by the variation 
of the location and form of the aerodynamic panel forces. For the strip theory 
representation, forces and moments are determined at a point on the elastic axis and at 
the midspan point of each strip. The PLEXSTAB, Doublet Lattice, and vortex 
spline-kernel function methods define forces that are constant over each panel and act 
at the panel centroids. The vortex spline-kernel function also defines moments at the 
centroid of each panel. The Rho unsteady kernel function method used a simplified 
integration of pressures to define lifts and moments on streamwise aerodynamic panels 
at the quarter chord. 


4.4 SOLUTION TECHNIQUES 

The method of solution of the equations of motion are identical whether the equations 
have either frequency dependent or constp.nt coefficients. The equations are solved for 
the generalized coordinates at discrete frequencies producing frequency response 
functions in terms of magnitude and phase angles for sinusoidal excitation functions. 


The excitation functions include aerodynamic forces resulting from either a continuous 
sinusoidal gust encounter (ref. 9) or a sinusoidal ailerbn oscillation. Once the 
generalized coordinate frequency responses were obtained, then the load equation 
frequency response functions were obtained by direct multiplication of the load 
coefficients by the generalized coordinate responses. Appendix B contains a more 
detailed discussion of the equations of motion and kad equation formulation. 


In utilizing the continuous sinusoidal gust excitation function, a comparison of the load 
parameters obtained by statistical methods produced a means of examining the 
adequacy of the aerodynamic methods for dynamic gust design loads with controls fixed. 
The two parameters of interest are A and Nq (refs. 11 and 12). 


A = 


. 1 1/2 

j lT(if2)l-cl>(S2)dn! 




(4) 


is the ratio of the root mean square (rms) value of the load to the rms value of the gust 
velocity. 


^ 2irA 


D 


11/2 


( 5 ) 


is the average number of times per unit distance that load crosses the value zero with 
positive slope; 


where: 


T(iii) = load frequency response function 
$(U) = gust spectrum 
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Note: The integrations in these expressions, though defined from zero to positive 
infinity, were actually performed from zero to 16 cycles per second (cps). 


The accepted gust spectrum for commercial airplane design purposes has been the 
Von Karman (fig. 8) which has the expression 


L 0-339 

[l +(1.33^) 


( 6 ) 


where; 

L = scale of turbulence, ft 
n = oi/v, spatial frequency, rad/ft 
Cl) = frequency, rad/s 
V = velocity, ft/s 

A and Nq can give an indication of the adequacy of the different aerodynamic- methods 
when dealing with continuous turbulance. However, these parameters are not able to 
give insight into the differences* in magnitude and phase angle prediction between the 
different aerodynamic theories. These differences can become critical in the design of 
active control systems for modal suppression, stability augmentation, or load 
alleviation. A clearer picture of the differences arising from the different methods can 
be obtained by examining the generalized coordinate and load frequency response 
functions resulting from gust and aileron excitation. 

Although two separate computer programs’*' were used depending on whether the 
coefficients:, to the equations were constant or nonconstant, the solution routines and 
output of each computer program are identical-only the ijiput routines differ. Thus, 
commonality was maintained wherever possible so that any differences were a result of 
only the aerodynamic methods. 


‘'‘Clemmons, H. E., A Power Spectral Digital Computer Program to Determine Dynamic Loads 
Due to Random Gusts-PSDSYS (TEV166)-User’s Guide, BCS-60235-1, June 1973. 

Sidewell, K. W., Solution. Program System (TEV114I for the Analysis of the Response of a 
Flexible Airplane to Vertical Turbulence, D6-29669TN, Vol. 1, March 1969. 



5.0 RESULTS AND DISCUSSION 


As discussed earlier, the aerodynamic methods used to generate equations of motion and 
load equations and which were compared are the: 

1 . Quasi-steady^‘ formulation of the FLEXSTAB steady state aerodynamics 

2. Quasi-steady formulation of the FLEXSTAB steady state aerodynamics with the 
instantaneous angle of incidence modified with Ktissner and Wagner indicia! lift 
growth functions 

3. Quasi-steady formulation of the FLEXSTAB steady state and low frequency 
approximation aerodynamics used both for response (motion induced) and gust 
aerodynamic forces 

4. Quasi-steady formulation of the FLEXSTAB steady state aerodynamics used for 
gust aerodynamic forces, and both FLEXSTAB steady state and low frequency 
approximation aerodynamics used for response (motion induced) aerodynamic 
forces 

Note: This method is the aerodynamics developed in FLEXSTAB and used in 
method 3 but with only the low frequency approximation aerodynamics on the 
response forces included. With the aerodynamic gust forces due to the flow 
incidence time rate of change eliminated, this method is theortically incomplete. 
However, it is useful for some stability and control analyses when evaluating the 
stability of the system. Thus, it was included as an aerodynamic method for use in 
comparison. 

5. Quasi-steady formulation of the strip theory aerodynamics corrected to match wind 
timnel steady state data, and the instantemeous angle of incidence modified with 
Kiissner and Wagner indicial lift growth functions 

6. Rho unsteady kernel function aerodynamics 

7. Vortex spline-kernel function unsteady aerod 3 mamics 

8. Doublet Lattice unsteady aerodynamics 

The solution of second-order differential equations becomes cumbersome and 
time-consuming when the coefficients of the equation are nonconstant with respect to 
frequency compared to an equation having constant coefficients. Unfortunately, the 
more exact unsteady aerodynamic methods listed as items 6, 7, and 8 produce frequency 
dependent equations (nonconstant coefficients) when used in dynamic loads analyses. 
Consequently, these methods have not been used extensively for dynamic loads 
predictions. 


*The quasi-steady formulation neglects the influence of the wake vorticies on the flow. This is 
equivalent to setting the Theordorsen function, C(k), equal to one (ref. 10). 


The accepted and most commonly used method has been the method listed as item 5. 
This method accounts for the unsteady aerodynamics in an approximate manner with 
lift growth functions but still retains its ease and speed in the equation solutions; the 
results are generally conservative for use in the design process. In addition, the 
theoretical aerodynamics can be modified in a simple manner with wind tunnel or flight 
test results. Unfortunately, this method is not suitable for aimlanes with low AR lifting 
surfaces. 

An analogous method to the strip theory (method 5) is the FLEXSTAB panel-type 
steady state aerodynamics modified with Kussner and Wagner indicial lift growth 
functions (method 2) to approximate unsteady aerodynamics. This method is applicable 
to airplanes with low AR lifting surfaces but has a disadvantage that it cannot 
presently be correctly modified to reflect wind tunnel or flight test results. 

The aerodynamic methods listed as items 1, 3, and 4 were considered in order to 
conclusively determine their adequacy or inadequacy for use in any type of dynamic 
loads analyses. The t 3 rpe of dynamic analysis that was performed with each 
aerodynamic method is shown in table 1. Note, for the oscillating aileron excitation, 
methods 3 and 4 are equivalent since these methods only differ in their approximation 
of the unsteadiness caused by the gust-forcing function. Method 1 was not run with the 
oscillatory aileron excitation because of the bad agreement with the unsteady Rho 
method in the gust excitation case (see fig. 9), 

The relative accuracy of these aerodynamic methods is ascertained by using the Rho 
unsteady kernel function aerodynamics as a basis for comparison. The roots of the 
characteristic equations for each of the aerodynamic theories, which yield constant 
coefficients, are shown in table 2. Since the other aerodynamic theories are frequency 
dependent, the roots of the equation vary with frequency and cannot be compared 
directly. 

Certain points of interest arise in studying the results shown in table 2. First, methods 
3 and 4 yield the same results. This should be expected since the difference between the 
two methods lies only in the approximation of the unsteadiness produced by the gust 
force which does not enter into the problem when rooting the equations of motion. 
Second, although the predicted frequencies are all approximately the same, there is 
considerable variation in the predicted damping. Third, when comparing each method 
over the range of modes, clear trends do not appear; therefore, no definite conclusions 
can be determined as to the relative accuracy of the prediction of loads by examining 
roots alone. It becomes necessary to turn to the comparisons of power spectral density 
(PSD) load parameters, frequency response function magnitudes and phase angles 
predicted by each method. 

5.1 PSD DYNAMIC GUST LOADS ANALYSIS, CONTROLS FIXED 

Table 3 shows the PSD load parameters for gust analyses with control surfaces fixed, 
obtained using random harmonic analysis techniques. The FLEXSTAB steady state and 
low frequency appwximation aerodynamics (method 3) used for both response and gust 
aerodynamic forces produce PSD load parameters whose values diverge with increasing 
. frequency. This is primarily caused by the effect of the low frequency approximation on 


gust angle, ag. This produces a generalized coordinate acceleration frequency response 
function that increases with increasing frequency* Since load is a function of inertia 
forces which are in turn a function of coordinate acceleration, the divergence of the load 
frequency response functions and A and Nq can be expected* Figure 9, which has plots 
of the magnitudes of wing root bending moment and wingtip acceleration frequency 
response functions due to continuous sinusoidal gust, demonstrates the effect of 
exceeding the reduced frequency restrictions for the low frequency approximations 
which were previously stated in section 4,2.1. In this section, the frequency restrictions 
were stated as: 


k« 1 

( 7 ) 

klinki « 2/ARyi-M- 


For the present analysis case, the Mach number was 0.8 and the aspect ratio was 6.29. 
This results in the reduced frequency being restricted to a value much less than 0.56 
radians and which corresponds to a cyclic frequency of 2.7 cps. Thus, the consideration 
of frequencies above 2 cps becomes meaningless for this method and, therefore, its use 
was discontinued. 

Using the Eho unsteady aerodynamics as a basis for comparison, the rms load/rms gust 
(A) percent difference for each aerodynamic method and each type of load is shown in 
table 4. 

Typical comparison of load magnitude frequency response plots due to gust excitation 
for the various aerodynamic methods, except method 3, with Rho unsteady 
aerodynaniics is shown in figures 10a through 12a.* The corresponding comparison of 
the phase angle relationships is shown in figures 10b through 12b.** For completeness, 
all of the comparisons of frequency response plots (magnitude and phase angle! for the 
generalized coordinates and loads are shown in appendix C (figs. C-1 through C-18). The 
corresponding PSD comparisons of figures 10 through 12 are shown in figures 13 
tbrovr': 16 for the response due to the Von Karman gust spectrum. These figures are 
determined by the expression: 

‘^'oLit=IT|-^>(f2) (8) 


*See table 5 for a guide of plot contents. 

‘rdinates of the phase angle figures represent the difference in phase angles predicted by 
metfiod i minus the phase angle predicted by Eho unsteady aerodynamics (method 6), 
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where: 


|T| = magnitude of load frequency response functions (figs. 10 through 12) 

= Von Karman’s gust spectrum 
^^out “ spectrum of the response quantity 

The PSD plots for such gust input are weighted at the low frequencies because of the 
gust spectrum. 

From the values of the PSD load parameters (A and Nq) given in table 3 (obtained from 
integration of the PSD plots), it can be concluded that for PSD gust loads calculated for 
the controls fixed configuration only the Doublet Lattice unsteady aerodynamics 
(method 8), the vortex spline-kernel function unsteady aerodynamics (method 7), strip 
theory aerodynamics modified with Kussner and Wagner indicial functions (method 6), 
and the PLEXSTAB steady state aerodynamics modified with Kussner and Wagner 
indicial functions (method 2) produce fairly satisfactory results compared to the Rho 
unsteady aerodynamics. However, this is only because of the large attenuation cf the 
gust power at &e higher frequencies (fig. 8). The other aerodynamic methods produce 
differences in A which, in some cases, become too large for design purposes. This 
becomes apparent when looking at the load frequency response plots (figs. 10 
through 12); the comparison of each aerodynamic method with Rho unsteady 
aerodynamics varies from fair to good for the frequency range from zero to 
approximately 2 cps. 

From 2 cps and higher, the comparison of all the aerodynamic methods with Rho 
unsteady aerodynamics (based on the magnitude of the frequency response function) is 
poor except for the vortex splme-kernel function and the Doublet Lattice aerodynamic 
representations which compare very favorably with Rho unsteady aerodynamics. The 
load frequency response phase angle plots (figs. 10b through 12b) also show appreciable 
differences over the entire frequency range between all of the aerodynamic methods and 
Rho unsteady aerodynamics, except for the vortex spline-kernel function and Doublet 
Lattice aerodynamics which generally compare favorably. 

At the several frequencies where the phase angle comparisons of the vortex 
spline-kernel function and the Doublet Lattice aerodynamic methods are poor, the 
magnitude of the load responses is near zero; e.g., compare figure 10a with 10b. This 
discrepancy in the phase angles occurs because small differences in the magnitude of 
the near zero load responses can result in large phase angle differences. However, since 
the load response magnitudes are small at these frequencies, it is expected that they 
would have little effect on an active control analysis regardless of the phase angle 
differences. 

5.2 DYNAMIC LOADS ANALYSIS, OSCILLATING CONTROL SURFACE 

Typical comparison of load magnitude frequency response plots due to an oscillating 
inboard aileron for aerodynamic methods 2, 4, 7, and 8 with Rho unsteady aerodynamics 


(method 6) is shown in figures 16a through 18a/^ The corresponding comparison of the 
phase angle relationship is shewn in figures 16b through 18b, For completeness, all of 
the comparisons of frequency response plots {magnitude and phase angle) for the 
generalized coordinates and loads are shown in appendix C (figs, C-19 through C-36). 

Prom the load frequency response plots (both magnitude and phase angle), only the 
results from the vortex spline-kernel function method and the Doublet Lattice method 
agree satisfactorily with the Rho unsteady aerodynamic results. Even with these 
unsteady aerodynamic methods, the agreement varies from excellent to only acceptable. 
It is apparent from these figures and those from section 5.1 that the aerodynamics are 
much more sensitive to an oscillating control surface than to a gust. 

Figures 16 through 18 also show that the agreement of the FLEXSTAB results 
(methods 2 and 4) with the Rho unsteady results is very poor when using an oscillation 
control surface. Examination of these figures shows a number of interesting occurrences 
concerning the three unsteady aerodynamic methods (methods 6, 7, and 8). First for 
oscillatory aileron excitation, the Doublet Lattice method, when compared to the Rho 
method, predicts higher frequency response function magnitudes for frequencies less 
than 3 cps; but for frequencies above 3 cps, the Doublet Lattice method underestimates 
the magnitudes. A study conduced at Boeing^''^ showed that the underestimation was a 
result of the paneling choosen. 

The Doublet Lattice method is a panel-type aerod 3 mamic method and is sensitive to the 
paneling scheme used. It was found in the Boeing study that the number of chordwise 
divisions must increase as the upper reduced frquency of interest increases in order to 
minimize potential error. For example, if the highest reduced frequency of interest were 
k = 0,5, then four chordwise divisions would be sufficient; whereas, if the highest k 
value of interest were k - 1.0, then eight chordwise divisions would be required* As 
shown in figure 5, six chordwise divisions were used over the majority of the wing in 
this study. Using the study of reference 4 as a guide, the highest k value which could be 
analyzed without considerable error is 0.7. In terms of the flight condilion used in this 
study, a k value of 0,7 is equivalent to a cyclic frequency of 3.4 cps. Therefore, the 
disagreement above 3 cps can be expected and its cause shown to be a result of the 
paneling scheme used. 

An examination of figures 10 through 12^ which are the frequency response functions 
for continuous gust excitation, shows that this discrepancy between Rho and Doublet 
Lattice above 3 cps does occur but is not as marked. This too is reasonable since a 
continuous gust excitation will tend to excite the higher modes less than a concentrated 
force such as an oscillating control surface. 


oscillatory control excitation, methods 3 and 4 are equivalent (see discussion in sec. 5.0) 
Method o was omitted because of the difficulty in modeling the oscillatory control using strip 
theory aerodynamics. Method 1 was left out because of its lack of agreement with Rho 
aerodynamics when using gust excitation. 

‘'"■Richard, M. and Kroll, R. L, ‘Ximitations of the Doublet Lattice Me thod--Appli cation to a High 
Aspect Ratio Wing/* Boeing document D6-42395. 1975. 


The reason for the low frequency disagreement between Rho and Doublet Lattice when 
using an oscillatoiy control ej:citation can only be a matter of conjecture at this time. 
Again, an examination cf figures 10 through 12 shows that Rho and Doublet Lattice 
agree extremely well in the low frequency region when using a gust excitation; yet for 
an oscillatory cohtrtil excitation (figs. 16 through 18), the agreement is not as good. 
Control surface deflections introduce high pressure gradients in the region of control 
surface boundaries. To model these gradients sufficiently with panel-type aerodynamics 
requires a fine grid in the control surface boundary region. The panel scheme for the 
osciUatory control surface case, shown in figure 5 , may not have been changed as much 
as required in the region of the control surface. 

A second point of interest in examining figures 16 through 18 is that the vortex 
spline-kernel function aerodynamics (method 7) consistently overestimates the response 
magnitudes. Looking to figures 10 through 12 for the gust excitation case shows that 
the agreement between Rho and this method is very good over the entire frequency 
range. Because the paneling scheme for this method (shown in fig. 2) was not changed 
between the two cases, it seems safe to assume that much of the disagreement is caused 
by the type of paneling scheme used. This assumption is supported by reference 2, 

The third point of interest is the phase angle differences shown in figures 16b through 
18b. In general, the phase angles^ agree quite well over the entire frequency range. The 
large sharp peaks and abrupt sign changes in phase difference curves are not 
significant, since they occur at very small magnitudes and are associated with the 
problem of division by small numbers. 

A final note of interest concerns the amount of aerodynamic damping predicted by the 
various methods when compared to Rho. Prom examination of the response plots due to 
gust and, even more so, due to oscillatory aileron, all aerod 3 mamic methods above 3 cps 
in general predict less aerodynamic damping than Rho as seen by their overestimation 
of response magnitudes. The approximate unsteady aerodynamic methods tend to 
predict less damping than the unsteady aerodynamic methods. 

A complete summary of figures 9 through 18 is given in table 5. 


6.0 CONCLUSIONS AND RSCOMMJiijNDATIONS 


Applicability of the various aerodynamic methods for use in dynamic loads analyses was 
ascertained by comparing results to those obtained using the Kho unsteady kernel 
function aerodynamic method which was used as a basis for comparison. The sensitiviiy 
of the results to the use of some of the various aerodynamic methods (methods 2, 5, 6, 
and 7) is not as critical when calculating controls fixed PSD gust loads versus loads due 
to active controls. This is because of the large attenuation of the gust-forcing function 
power at the higher frequencies. However, even with this relaxed sensitivity to the 
aerodynamic representation, only the PSD gust loads (A) calculated using the vortex 
spline-kernel function (method?) and the Doublet Lattice aerodynamics (methods) 
differed less than 10%* from those calculated using Rho unsteady aerodynamics 
(method 6). RMS loads (A) calculated with strip theory and FLEXSTAB steady state 
aerodynamics modified with Kiissner and Wagner functions differed from -12% to +18%. 
The remaining aerodynamic methods, FLEXSTAB steady state (method 1), FLEXSTAB 
steady state plus low frequency approximation on response aerodynamics only 
(method 4), and FLEXSTAB steady state plus low frequency approximation on both the 
gust and response aerodynamics (method 3), produce A’s which varied much above 20% 
with results from the latter method actually diverging. 

When active controls are used to influence elastic structural motions, none of the 
quasi-steady state aerodynamic representations are satisfactory even with unsteady 
effects represented with indicial lift growth functions. To properly represent the 
aerodynamics resulting from an oscillating control surface, more exact unsteady 
aerodynamic methods must be employed. 

There are some advantages in using tlie approximate methods in representing unsteady 
aerodynamics. These methods arc- less expensive to use than the fiill unsteady methods 
because they yield equations of motion which are easier to formulate and solve. In 
addition, the modal responses tend to be less damped than those using unsteady 
aerodynamics and result in more conservative loads. The disadvantages are that the 
response and loads, magnitudes and phase angles are inaccurate at higher frequencies 
and cannot be used for active control design with any amount of confidence. 

This study used an advanced state-of-the-art unsteady aerodynairiic method as a basis to 
determine the validity of various aerodynamic theories. It is assumed that this theory 
satisfactorily matches the actual aerodjmamics on a planar surface. There have been 
very few oscillatory aerodynamic wind tunnel tests performed on lifting surfaces to 
credit or discredit any exact unsteady aerodynamic methods at oscillatory frequencies 
other than zero. ' 


^'The vortex spline-kernel function method results differ by «3% from the Rho aerodynamic 
method. 


A logical extension of this present study would be to perforin a theoretical dynamic 
vertical gust analysis on an airplane employing total airplane flexibility and a 
sufficient number of modes to adequately cover the gust bandwidth as deffiied from 
flight test results. The theoretical gust analysis load frequency response functions 
would then be compared to the load frequency responses obtained from gust loads 
survey flight test results. Several aerodynamic theories could be used in the theoretical 
vertical gust analysis to verify their adequacy in gust loads predictions on the total 
airplEine and indications of their adequacy for active controls representations. The 
aerodynamic theories are (1) Doublet Lattice unsteady aerodynamic method and 
(2) quasi-steady formulation of the FLEXSTAB steady state aerodynamics modified with 
Kiissner and Wagner indicial lift coefficients. 


Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle, Washington 98124, July 1976 


Table 1,— Type of Dynamic Analysis Performed With Each Aerodynamic Method 




Gust excitation function 

Oscillating aileron 
•• excitation function 

Method 

Aerodynamics based on 

Loads and responses 

Loads and responses 



Magnitude 

Phase angle 

Magnitude 

Phase angle 

1 

FLEXSTAB steady state • 

- / 

y ^ 



2 

FLEXSTAB steady state rnodified with Kiissner and Wagner Functions 

y 


/ 

y 

3^ 

FLEXST AB steady state plus low frequency approximation 
used both for response and gust forces 

V 

“ • 



4" 

FLEXSTAB steady state for gust forces and both FLEXSTAB ’ _ 

steady state plirs low frequency approximation for response forces 

y ■ 

r 


y 

5 

Strip theory corrected to wind tunnel daca and modified 
with Kussner and Wagner functions 

y , 

y 

\ 



6 

Rho unsteady kernel function 

y 

y 

y 

y 

7. 

Vortex spline-kernel function 

y 

y 

y 

y 

3 

Doublet lattice j 

V . 

- y 

y . 

y 


®The differencG between methods 3 and 4 fs only how the gust excitation force is used. 
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raib/e2.— jE^7uaf/b/7Sof /Wof/o/7 Ftoots^ 


Aerodynamic 

method 

Modes 

Basic 
structural 
frequencies 
and damping 
(in vacuo) 

FLEXSTAB aerodynamics 

5 

1 

2 

3 

4 

Strip theory 
corrected to wind 
tunnel data with 
Kussner and Wagner 
k = 0 

Steady state^ 
without 

Kussner and Wagner 

Steady state 
with 

Kussner and Wagner 

Steady state + low 
freq apiDrox aero 
on response and 
gust aero 

Steady state low 
freq approx aero 
on response aero 

fn.cps 

f = C/Cg 

f„, CPS 

f = c/Ct, 

fn-cps 

f = c/Cg 

fn-cps 

? = c/Cp. 

fn. ■=Ps 


fp, cps 

f •-= C/Cg 

Translation 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Pitch 

0.0 

0.0 

0.233 

0.479 

0.235 

0.475 

0.234 

0,489 

0.234 

0.489 

0.261 

0.46T 

First elastic mode 

1.00 

0.0 

1.16 

0.097 

1.17 

0.085 

1.17 

0.073 

1.17 

0.078 

1.16 

0.093 

Seuoiid elastic mode 

2;i5 

0,0 

2.49 

0,061 

2.48 

0.033 

2.51 

O.QGO 

2.51 

0.060 

2.35 

0.058 

/.Third elastic mode 

3.04 

0.0 

3.11 

0.059 

3.13 

0.045 

3.07 

0.082 

3.07 

0.082 

3.08 

0.056 

Fourth elastic mode 

4.16 

0.0 

4.13 

0.007 

4-14 

0.007 

4.09 

0.024 

4,09 

0.024 

4.12 

0.008 

Fifth elastic mode 

6.17 

0.0 

6.28 

0.028 

6.28 

0.017 

6.26 

0.036 

6.26 

0.036 

6.13 

0.032 

Sixth elastic mode ^ 

8.25 

P.O 

8.30 

0.018 1 

8.31 

0.012 

8.20 

0.035 

8.20 

0.035 1 

8.18 

0.026 


^Baseci on large subsonic jet airplane; weigh! = 753 000 lb; altitude = 20 QOD ft; IVIach = 0.8; eight free-free mode (\/ertical translation, pitch, and six elastic wing modes). 
^Quasi'Steady approach (see sec. 5, OJ: 



Table 3.— PSD Load Parameter Comparison^ 



FLEXSTAB aerodynamics 

Aerodynamic 

method 

1 

2 


4 

Steady state 
without 

Kussner and Wagner 

Steady state 
with 

Kussner and Wagner 

Steady state + low 
freq approx aero 
on response and 
gust aero 

Steady state + low 
freq approx aero 
on response aero 

A 

Wo 


Wo 

A 

Wo 


Wo 

Shear, lb 

Load 

per 

fps 

Zero 
crossings 
per feet 

Load 

per 

fps 

Zero 
crossings 
per feet 

Load 

per 

fps 

Zero 
crossings 
per feet 

Load 

per 

fps 

Zero 
crossings 
per feet 

2072 

0.00317 

1924 

0.00258 

5002 

0.01144 

1943 

0.00237 

Root 

Outbd of inbd nac 

1568 

0.00146 

1599 

0.00138 

2174 

0.01033 

1705 

0.00139 

Outbd of outbd nac 

551 

0.^’0166 

543 

d00135 

1246 

0.01094 

576 

0.00157 

Bending moment, 
in-lb X 10® 

Root 

1.1100 

0.00196 

1.1020 

0.00163 

2.629 

0.01111 

1.1650 

0.00147 

Outbd of inbd nac 


0.00137 

0.5620 

0,00128 

1.229 

0.01062 

0.6090 

0.00131 

Outbd of outbd nac | 


0.00280 1 

0.0837 

0.00230 

0.197 

0.01103 

0.0926 

0.00234 

Acceleration, g s 
Wingtip 

0.1870 

0.00690 

0.1480 

0.00580 

0.2650 

■ 

0.00912 

0.1720 

0.00694 

Outbd of outbd nac 

0.0519 

0.00417 

0.0505 

0.00328 

0.0676 

0.00652 


0.00371 

Outbd of inbd nac 

0,0232 

0,00603 

0.0198 

0.00478 


0.00851 

0.0213 

0.00507 

Root 

0.0138 

0.00268 , 

0.0131 

d00208 


0.00316 1 

0.0134 

0.00210 


^Based on large subsonic jet airplane: weight - 753 000 lb; altitude = 20 000 ft; Mach = 0.8; eight free-free mode 
(vertical translation, pitch, and six elastic wing modes). 

“Loads have not converged. 

Note: cutoff frequency is 16 cps. 
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Table 3.— (Concluded) 


Aerodynamic 

method 

Loads 

5 

6 

7 

8 

Strip theory 
corrected to wind 
tunnel data with 
ICussner and Wagner 
k = 0 

Rho 

unsteady 

aerodynamics 

Vortex spline- 
kernel function 
nonsteady 
aerodynamics 

Doublet lattice 
unsteady 
aerodynamics 

Shear 

Root 

A 

^>0 

A 

No 

A 

No 

A 

No • 

Load 

per 

fps 

Zero 
crossings 
per feet 


Zero 
crossings 
per feet 

Load 

per 

fps 

Zero 
crossings 
per feet 

Load 

per 

fps 

Zero 
crossings 
per feet 

1830 

0.00259 

1708 

0.00226 

1705 

0.00226 

1686 

0.00221 

Outbd pf inbd nac 

1474 

0.00127 

1513 

0,00116 

1481 

0.00116 

1558 

0.00114 

butbd of outbd nac 

B14 

0.00119 

516" 

0.00112 

515 

0.00114 

536 

0,00105 

Bending moment, 
in-lb X 10 ® 

Hoot 

0.9350 

0.0020 

1.05 

0.00132 

1.03 

. 

0*00134 

1.07 

0,00135 

Outbd of inbd nac i 

0.5270 

0.0013 

0.547 

0.00125 

0.543 

0.00122 

0.567 

0.00120 

Outbd of outbd nac 

0.0841 

0.0018 

0.075 

0.00225 

0.077 

0.00215 i 

0.081 

0.Q0195 

Acceleration, g s 
Wingtlp ' 

0,1170 

0.00507 ! 


0.00687 

0.1316 

0.00676 1 

0.1294 

O.OD653 

Outbd of outbd nac i 

a 0429 

0.00266 

HEESS9 

0.p03“51 , 

0,0431 

0.00346 i 

0.0447 

0.00335 

Outbd of inbd nac ' 

0.0177 

0.00409 

0.0179 

0.00673 1 

0.0175 

. 0,00656 

0.0175 

0.00631 

Root 

0,0136 

0.00198 ! 

0.0126 

0.00188 

0.0122 

0.00195 

6.0123 



Note: cutoff frequency is 16 cps 







































































































Table 4.— Percent Difference in Load of Each Aerodynamic Method Compared to Rho Unsteady Aerodynamic Method 


Method 

ABradynamic method 

A percent difference from Rho {method 6) 

Wing shear at 

Wing bending moment at 

Wing acceleration at 

Root 

Inbd 

nacelle 


Root 



Tip 


Inbd 

nacelle 

r- 

Root 

1 

FL EXSTAB steady state 

2T0 

4.0 

7.0 

6.0 

1.0 

23.0 

40.0 

19.0 

30.0 

10.0 

2 

FLEXSTAB steady state with Ktissner and 
Wagner incliciat lift growth functions 

13,0 

6.0 

6.0 

5.0 

3.0 

18.0 

11.0 

16.0 

11.0 

4.0 

3 

FLEXSTAB steady state plus low frequency 
approximation on response and gust aerodynamics 

193,0 

116.0 

141.0 

150.0 

125.0 

163.0 

S8.0 

55.0 

61.0 

13.0 

4 

FLEXSTAB steady state and low frequency 
approximation on response aerodynamic 

14.0 

13.0 

12.0 

11.0 

11.0 

23.0 

29.0 

30.0 

19.0 

6.0 

5 

Strip theory with Kiissner and Wagner i 

indicial lift growth functions 

7.0 

-3.0 

0.0 

-11.0 

-4.0 

12.0 

-12.0 

-2.0 

-1.0 

8.0 . 

. 7 

Vortex splipe-kernel function 

0.0 ' 

-2.0 

0.0 

■HBl 

-1.0 

3.0 

-1.0 i 

-1.0 

-2.0 

-3,0 

8 

Doublet lattibe 

-1.0 

3.0 

4.0 

2.0 

4.0 

8.0 

-3.0 

2.0 

-2.0 

-2.0 


► 



























































































Table 5,— Contents of Companson Plots tn Figures 9 Through 18 
Response Description (Two Rigid Body Modes and Six Elastic Modes) 


Response variable 


Frequency response due 
to gust ejccitation 


ft/lethods 
3 and 6 


All methods 
excepts 


Magnitude Magnitude 


>/ 


Power spectral 
density due to 
Von Karnnan Input 


All methods except 3 


7 

v/ 

V 

V 


V 


Frequency response due 
to oscillating inboard 
aileron 


Methods 2, 4, 6, 7, and 8 


Magnitude 
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FREQUENCY. CYCLES PER SECOND 


FREOUCNCY. CYCLES PER SECOND 


Figure 9.— Load Frequency Response Function Magnitudes Due to Gust Excitation 
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Figure 10.— Comparison of Wing Root Shear Frequency Response Function Due to Gust Excitation 
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Figure 10,— (Concluded) 
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Btuiso ••>5 s-s' ; '’es»:'.s£. : 





(a) Magnitudes 

Figure 1 1 .—Comparison of Wingfioot Bending Moment Frequency Response Function 
Due to Gust Excitation 
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Phase angle relative to rho AEROovNAr^ics. RAO. Phase angle relative to r-2 aercz^nah: 





(b) Phase Angles 

Figure 11.— (Concluded) 
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Figure 12.-Comparison of Wingtip Acceleration Frequency Response Function Due to 
Gust Excitation 
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Figure 12.— (Concluded) 
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Figure 13.— Wing Root Shear Power Spectra! Density 
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Figure 16.— Comparison of Wing Root Shear Frequency Response Function 
Due to Aileron Excitation 
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(a) Magnitudes 




(b) Phase Angles 

Figure 17.— Comparison of Wing Root Bending Moment Frequency Response Function 
Due to Aileron Excitation 
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Figure 18.— Comparison of Wingtip Acceleration Frequency Response Function Due to 
Aileron Excitation 



APPENDIX A 


STRUCTURAL MODEL 


In order to simplify the structural model but maintain a meaningful dynamic analysis, 
the airplane structural model was represented by only the wing elastic degrees of 
freedom. The body, empennage, and nacelles were considered to be rigid. The wing was 
modeled as a finite number of lumped masses, rotary inertias, and beam segments 
(fig. A-1). The lumped masses were considered to be attached to the elastic axis with 
rigid links perpendicular to the elastic axis. The flexibility matrix for this model was 
obtained and combined with the mass representation to form an eigenvalue problem. 
This problem was solved for free-free mode shapes and frequencies. The generalized 
mass and stiffness matrices were obtained for the Tirst eight free-free modes which 
consisted of two rigid body modes (vertical translation and pitch) and six wing elastic 
modes. This representation of the generalized mass and stiffness was used for each 
aerodynamic method in formulating the equations of motion. 

The elastic mode shapes obtained from the solution of the eigenvalue problem consisted 
of vertical bending displacement ((f>BdR^ torsion (</>Tf>rsl bending slopes at each 
point on the elastic axis where the ma.s.s was considered to be attached by rigid links 
(figs. A-2 through A-4).* 

These mode shapes were then interpolated to the necessary aerodynamic reference 
points required by each aerodynamic method. Figures A-5 through A-11 show the grid 
used for the !7LEXSTAB aerodynamics and the interpolated mode shapes at the 
associated aerodynamic reference points. The interpolation routine employed a cubic 
spline in arc lengths along a continuous planar curve (motion axis)** with mapping 
from the motion axis along defined reference lines which are assumed rigid with respect 
to displacement and rotations at the motion axis. In this case, tht' motion axis and 
elastic axis are the same. This interpolation routine was used to interpolate the mode 
shapes from the structural node points to the associated aerodynamic reference points of 
each aerodynamic method. 


*ln figures A-2 through A-4. the ordinate for torsional displacement is read as ffrxors ^ so 
that for a value of -0.4 read from the plots, the actual d>Tors ~ -0.4/100. 

**A motion axis is a line along which the displacement due to rigid airplane and elastic motion is 
defined. 





NORMALIZED MADE SHAPE NORMALIZED MADE SHAPE 

^Bdg'*^Tors^’00 ^Bdg'^Tors ^ 



First Elastic Mode 
Freq » 1.00 cps 


ELASTIC AXIS WING SPAN STATION, INCHES' 


Elastic axis wing span is the cummulative distance along 
the elastic axis starti^^g at the wing root. 

*AII mode shapes normalized by tip deflections are units of 
in./in. and^x^rc rad/in. 


Second Elastic Mode 
Freq = 2.15 cps 




ELASTIC AXIS WING SPAN STATION, INCHES 


Figure A-2.— First and Second Elastic Mode Shapes on the Elastic Axis 




NORMALIZED MADE SHAPE NORMALIZED MADE SHAPE 

**Bdg' "^Tors ^ ^ “^Bclg' ^Tors ^ 


Third Elastic Mode 
Freq = 3.04 cps 


100 2 300^ 4 _ 500 6 700 8/ 900 10 1100 


ELASTIC AXIS WING SPAN STATION, INCHES 


Fourth Elastic Mode 
Freq = 4.16 cps 


100-^ 300 4 500 6 700 v_8_^900 10 1100 ^Tors 


ELASTIC AXIS WING SPAN STATION, INCHES 


Figure A-3,— Third and Fourth Elastic Mode Shapes on the Elastic Axis 



APE 



Fifth Elastic Mode 
Freq = 6.17 cps 


ELASTIC AXIS WING SPAN STATION, INCHES 


Sixth Elastic Mode 
Freq = 8.25 cps 


ELASTIC AXIS WING SPAN STATION, INCHES 


Figure Fifth and Sixth Elastic Mode Shapes on the Elastic Axis 



NCHES 


Plotting 



Figure AS.— Modal Geometry Definition Used With FLEXSTAB Aerodynamic Grid 
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APPENDIX B 


EQUATIONS OF MOTION 
AND LOAD £ ;UATIONS FORMULATION 

The equations of motion were formulated by using the Lagrangean approach. Two basic 
forms were developed depending on the aerodynamic representation: (1) constant 
coefficient equations formed by utilizing stea V state aerodynamic methods and 
approximating unsteady aerodynamics with either Kiissner and Wagner indicial lift 
growth functions or the FLEXSTAB low frequency representation and (2) nonconstant 
coefficient equation with respect to frequency formed utilizing the more exact unsteady 
aerodynamic methods. 

The generalized structural inertia and stiffness coefficients were obtained from the 
structural representation defined in appendix A. Although these coefficients are 
constant with respect to frequency, they were used in conjunction with both of the basic 
forms defined above. 

The aerodynamic methods which produced constant coefficient equations for the 
generalized aerodynamic forces were:* 

1. Quasi-steady formulation of the FLEXSTAB steady state aerodynamics 

2. Quasi-steady formulation of the FLEXSTAB steady state aerodynamics and the 
instantaneous angle of incidence modified with Kiissner and Wagner indicial lift 
growth functions 

3. Quasi-steady formulation of the FLEXSTAB steady state and low frequency 
approximation aerodynamics used both for response and gust aerodynamic forces 

4. Quasi-steady formulation of the FLEXSTAB steady state aerodynamics used for 
gust aerodynamic forces, and both FLEXSTAB steady state and low frequency 
approximation aerodynamics used for response aerodynamic forces 

Note: This method is the aerodynamics developed in FLEXSTAB and used in 
method 3 but with only the low frequency approximation aerodynamics on the 
response forces included. With the aerodynamic gust forces due to the flow 
incidence time rate of change eliminated, this method is theoretically incomplete. 
However, it is useful for some stability and control analyses when evaluating the 
stability of the system. Thus, it was included as an aerodynamic method for use in 
comparison. 


*Note: These methods are numbered in the same order as in table 1 and as in section 5.0. 


5. Quasi-steady formulation of strip theory aerodynamics corrected to match wind 
tunnel steady state data, and the instantaneous angle of incidence modified with 
Kiissner and Wagner indicial lift growth functions 

The other basic form of equations, nonconstant coefficients formulated using the exact 
unsteady aerodynamic methods, consisted of the: 

6. Rho unsteady kernel function aerodynamics 

7. Vortex spline-kernel function unsteady aerodynamics 

8. Doublet Lattice unsteady aerodynamics 

The equations of motion* have the following form for the constant coefficient equations 
(refs. 8 and 13): 

lMll(q} + (M3U‘l}+ (M4l{q}*‘l>+ IM5I {?}*«!>= {CjlF*./' (B-1) 

A variation of equation (B-1) for the frequency dependent coefficients from the unsteady 
aerodynamic methods (refs. 14 and 15) is shown in equation (B-2). 

lM3l(q}+[M4 + M|l{q} + [M5l{q} = {C2}F (B-2) 

where: 


[Ml] 

[M3] 

[M4] 

[M5] 

[c.l 

{C2} 

F 

cl> 



generalized structural stiffness matrix 
generalized structural inertia matrix 
generalized aerodynamic stiffness matrix 
generalized aerodynamic damping matrix 
generalized forcing function coefficients (real) 
generalized forcing function coefficients (complex) 
forcing function 
Wagner indicial lift function 
Kussner indicial lift function 


*The equations are formulated in the inertial axis system for the initial condition of straight and 
level flight. 


indicial convolution 


q = generalized coordinates 

The equation of motion formulation for aerodynamic methods 1 through 4 used the "s” 
plane formulation as given in reference 7 with the exception of the low frequency 
approximation aerodynamics, which requires the rate of change of the downwash (W), 
given by: 

{w}= I0’l{q} +~l0l{q} (B-3) 


as well as the downwash (W) which is required when using the FLEXSTAB steady state 
aerodynamic force formulation (ref* 1). 

The strip theory aerodynamics (method 5) also used the **s’* plane formulation. The 
aerodynamic lift was formulated using two-dimensional incompressible flow. Steady 
state three-dimensional lift was then obtained by modifying the two-dimensional lift 
coefficients with compressibility factors (Fc), wind tunnel test data, and a static 
induction matrix based on the method developed in reference 6. From this, the 
generalized aerodynamic force coefficient matrices were formed. 

The vortex spline-kernel function aerodynamics produced an AIC matrix which relates 
panel downwash to panel pressure, which was used as in the method presented in 
reference 7. The Rho unsteady kernel function and the Doublet Lattice unsteady 
aerodynamic representations could not be used in this manner, as explained in the 
following paragraph* 

The output from the Rho kernel function and Doublet Lattice aerodynamic program 
consisted of complex matrix coefficients which yielded generalized forces and surface 
pressures when multiplied by the generalized coordinates (q)* 


iOl = (lORel+ilO|ml){q} 

{P}=(lPRel+i|P|„l){q) 


where: 


{Q}, {P} = complex generalized forces and surface pressures, respectively 

[QRe]t [Qlin]= real and imaginary parts of the generalized force coefficients obtained 
from the aerodynamic programs 

[pRe]» [Pim] = real and imaginary parts of the pressure coefficients obtained from the 
aerodynamic programs 


= generalized coordinate 


The solution program required that the generalized force coefficients be coefficients of q 
and q. Therefore, using the sinusoidal relationship 

q = icjq 

the generalized forces and surface pressures were written as: 

{Ol=IORel{<l)-H/<^IQ|„,l(ql 

|P|=IPReHq| + IMP|„Mq} 

Further modification was made to the generalized force and surface pressure coefficient 
matricies to account for nondimensionalizing. 

In addition, because of the computer time required to calculate the unsteady 
aerodynamic generalized forces and pressures at each frequency, only a limited number 
of base frequencies were calculated (20), and the matrices were then linearly 
interpol^*^ed to a larger number of frequencies (210) so as to adequately define the 
transfer functions once the equations were solved. 

The load equations have a form similar to the equations of motion: 

load= lM3l{q} + (M4l{q}*4> + [M5l{q}^cli + {c,}F* (B-7) 

A variation of equation (B-7) for the frequency dependent coefficients from the unsteady 
aerodynamic methods is shown in equation (B-8), 

load = [R3l{c7 }+(M 4 l{q} +|M5l{q}+{C2}F (B-8) 

where: 

M3 = inertia force coefficients 

M4 = aerodynamic stiffness force coefficients 

M5 = aerodynamic damping force coefficients 

Cl = real matrix coefficients of loads due to the forcing function 

C2 = complex matrix coefficients of loads due to the forcing function 

F = forcing function 

Although the computer program used to solve the constant coefficient equations differed 
from that used to solve the frequency dependent coefficient equations, there existed 
commonality betw^een the two programs in the solution and output routines. Only the 
input routines were different. The computer program used for the frequency dependent 
coefficient equations allowed input of matrices at various base frequencies and then 


(B-5) 


(B-6) 


linearly interpolated the matrices for a larger number of solution frequencies. In 
addition, complex matrices for the forcing functions could be used as input. 

The acceleration load equation matrix coefficients for all of the aerodynamic methods 
were formulated and generated as defined in reference 7. The shear and bending 
moment load equation matrix coefficients were also formulated and generated as 
defined in reference 7 for the aerodynamic methods 1 through 4, 7, and 8, and for the 
coefficients representing the inertia forces used with all of the aerodynamic methods. 
The strip theory (method 5) represented aerodynamic forces and moments at a point on 
the elastic axis and at the midspan point of each strip. These forces and moments were 
summed to calculate the shear and bending moment matrix coefficients at the specific 
load stations. 

In order to calculate the shear and bending moment matrix coefficients using the Rho 
unsteady kernel (method 6), it was necessary to perform a chordwise integration of the 
chordwise pressure distributions reported by Rho, since this aerodynamic method is a 
collocation point method and not a panel-type method. This procedure, using a 
trapezoidal integration scheme, produced lifts and moments at the quarter point of each 
chord where pressures were reported. These lifts and moments were then integrated to 
the required load stations to produce the appropriate shear and bending moment matrix 
coefficients. 

Parameters used in the equations of motion, load equations, and solution routines for 
the Kiissner and Wagner indicial lift growth functions and gust parameters were: 

Gust spectrum Von Karman 

Gust velocity 1 fps 

Scale of turbulenje 2500 ft 

Kiissner indicial lift growth coefficients (where applicable): 

V V V 

-0.58— t -1.45 — t -6.0 — t 

1.0 -0.488 e -0.272 c -0.193e (B-9) 

Wagner indicial lift growth coefficients (where applicable): 

V 

-0.762-^t 


0(t)= 1.0-0.361 e 


(B-10) 


APPENDIX C 


GENERALIZED COORDINATES 
AND LOAD FREQUENCY RESPONSE FUNCTIONS 

The generalized coordinates and load magnitude frequency response plots for the 
solutions of the equations identified in appendixes A and B are shown in figures C-la 
through C-36a. The corresponding phase angle difleren<’e (t^meUiodi * frequency 

response plots are shown in figures C*lb through C-36b. 

The solutions of the equations were separated into two types of analysis consisting of: 
(1) dynamic gust loads analysis with controls fixed and (2) dynamic loads analysis 
performed with the excitation function defined by an oscillatory control surface (inboard 
aileron). The transfer functions for the dynamic gust loads analysis with controls fixed 
are shown in figures C-la through C-18a for magnitudes and figures C-lb through C-18b 
for phase angles. The transfer functions for the dynamic loads analysis performed with 
an oscillatory control surface excitation function are shown in figures C-19a 
though C-36a for magnitudes and figures C-19b through C -36b for phase angles. A 
content summary of figures C-1 through C-36 is given in table C-1. 

The Rho unsteady kernel function aerodynamic method is presently the theoretically 
most exact unsteady aerodynamic representation. Consequently, it was used as a basis 
for comparing the accuracy of the other aerodynamic methods, and its transfer functions 
are presented in each figure for comparison purposes. 
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Table C- 1.— Content of Comparison Plots in Appendix C Figures 
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Figure C-W.— Comparison of Wing Shear Frequency Response Function at 
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Figure C-lh— Comparison of Wing Shear Frequency Response Function at 
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Figure C- 12.— Comparison of Wing Bending Moment Frequency Response Function 
at Wing Root Due to Gust Excitation 
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Figure C-12,— (Concluded) 
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Figure C- 13.— Comparison of Wing Bending Moment Frequency Response Function 
at Inboard Nacelle Due to Gust Excitation 
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Figure C- 15.— Comparison of Wing Acceleration Frequency Response Function at 
Wingtip Due to Gust Excitation 
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(b) Phase Angles 

Figure C- 15,— (Concluded) 
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Figure C-16. -Comparison of Wing Acceleration Frequency Response Function at 
Outboard Nacelle Due to Gust Excitation 
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Figure C- 16 , ’-'(Concluded) 




98 



(a) Magnitudes 

Figure C-1 7. —Comparison of Wing Acceleration Frequency Response Function at 
Inboard Nacelle Due to Gust Excitation 
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Figure C-17.- (Concluded} 
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(b) Phase Angles 

Figure C-18.-(Concluded) 
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Figure C-1 9.— Comparison of Gtneralized Coordinate Frequency Response Function 
Due to Oscillatory Aileron Excitation 


(a) Magnitudes 


X STT*:- .ran mxmti rw 

^L£r?TAf surr * ltw 


FREQUENCY. CYCLES PER SECOND 


STATE ■ITM Fra 


-U I 


FRECUENCY, CYCLES PER SECOND 


7 «l£FS7« S' 

X FLf*ST*c S' 
'?ra:j 


103 


□ V 



I I 



.• 4 (| 

i 


I I I 


0 5W0 

i > #0«TIl fUWCTlON 

□ JuuBtET lA’TlCI 


c»:.r:s ^ s- 


(a) Magnitudes 


1- 


‘ * \ 


*✓. - 
\ 


t 


j r. fiMJTi? ^TCA^ 
I A HEItTEF '•?!:» 




j "! ri 


I 

1 . 


• i... 


! 


•f-- 


FRECjJENCV« cycles per 5 EC 0:»0 


J ; • 


[... |... { 1 1 

- 



rv r :v: • 

. i . i ..... 

■ ‘ ‘ 

. - r«'. 

: j ' 


< I 

\ \ \ 

> wOLttLtT L-mr.E 


• * 1 

; J r : 

■ ! ; : ! i 

> 


: 1 - • • 




1 ' 'f 

1 ' : : ; ; i ! 


1 ' : i 

■ 





i I 


i.i'V*. i.*^»cs ctR 


(b) Phase Angles 

Figure C-20.— Comparison of Generalized Coordinate Frequency Response Function 
Due to Oscillatory Aileron Excitation 
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(b) Phase Angles 

Figure C-21 .-^Comparison of Generalized Coordinate Frequency Response Function 
Due to Oscillator/ Aileron Excitation 
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Figure C-22,— Comparison of Generalized Coordinate Frequency Response Function 
Due to Oscillatory Aileron Excitation 
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lb) Phase Angles 

Figure C-23.— Comparison of Generalized Coordinate Frequency Response Function 
Due to Oscillatory Aileron Excitation 
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(b) Phase Angles 

Figure C-24.— Comparison of Generalized Coordinate Frequency Response Function 
Due to Oscillatory Aileron Excitation 
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Figure C-25.— Comparison of Generalized Coordinate Frequency Response Function 
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Figure C-26,— Comparison of Generalized Coordinate Frequency Response Function 
Due to Oscillatory Aileron Excitation 
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(b) Phase Angles 


Figure C-27.— Comparison of Wing Shear Frequency Response Function at Wing Root 
Due to Oscillatory Aileron Excitation 
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Figure C-28.— Comparison of Wing Shear Frequency Response Function at Inboard 
Nacelle Due to Oscillatory Aileron Excitation 
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(b) Phase Angles 


Figure C-29.— Comparison of Wing Shear Frequency Response Function at Outboard 
Nacelle Due to Oscillatory Aileron Excitation 
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Figure C-30.— Comparison of Wing Bending Moment Frequency Response Function at 
Wing Root Due to Oscillatory Aileron Excitation 
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(b) Phase Angles 


Figure C-32.— Comparison of Wing Bending Moment Frequency Response Function at 
Outboard Nacelle Due to Oscillatory Aileron Excitation 
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(b) Phase Angles 

Figure C-33.— Comparison of Wing Acceleration Frequency Response Function at 
Wingtip Due to Oscillatory Aileron Excitation 


117 




PHfSE angle relative to RhO AEROOvNA!iICS» RAO« • ILiMK . RESPONSE.t . ,/10 OEG. AtL. T.OT 



O J*»SU- ■ nu'^. 

X VWTfl SPLl«*rf*?lE. F.iTTlWl 
□ j :.LT L*TTj:r 


V 




■A 


I 


/ 
/ ‘ 


,1 

S I 


A,/ 


V- 


(a) Magnitudes 




< , 


• - FiElSr&3 5Tra:» £TArr 

•; • Fa»jTi«! iTr:r:» T'r*f * 

APP*TK. 


•- - 

• -• 

» 

f « 


t j ’ * i 

lv.Vs/ 

/ 


<r t . ^ r— 

V A • 

• \ 

• 4 


t <. 

I t 

1 4 ' 1 ; 


* ’ * • i 

_i\ ‘ 


• ■ : r 

: . t { : 

t - : . .) 

’ ! \i 

,1 »! . 


i ; 1 r .1 

1 X 

1 ‘ • 

0 

... ! : i .K 

•» 1 . ; 

t t 1 4 ) t 

' ’ * i ' 

F • « 1 


UJ *' 

<r 

UJ 

< 

^ *• 





t 

1 

•— 1 

— 

t 

"■f 

V. 

.. , i» 

4 

It 


1 









I 

.::j 





; o ,rU:i£-CI!lJH£. Fi.«NCT:w> 

: ; i V Juoait ■ ia:t:u 


T' 


“I 

i 

1 




f i 


' KtH 




» ■ . i 

— - ! :{: 




FRECJcNCY. cycles PER SEC0N3 


FREOUENCY. CYCLES PER 5EC0N0 


(b) Phase Angles 

Figure C-34,— Comparison of Wing Acceleration Frequency Response Function at 
Outboard Nacelle Due to Oscillatory Aileron Excitation 
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Figure C-35.— Comparison of Wing Acceleration Frequency Response F unction at 
inboard Nacelle Due to Oscillatory Aileron Excitation 
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Figure C-36.—Coi nparison of Wing Acceleration Frequency Response Function at 
Wing Root Due to Oscillatory Aileron Excitation 
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